Effects of polarization on the transmission and localization of classical waves in weakly 

scattering metamaterials 
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We summarize the results of our comprehensive analytical and numerical studies of the effects 
of polarization on the Anderson localization of classical waves in one-dimensional random stacks. 
We consider homogeneous stacks composed entirely of normal materials or metamaterials, and also 
mixed stacks composed of alternating layers of a normal material and s metamaterial. We extend the 
theoretical study developed earlier for the case of normal incidence [A. A. Asatryan et al, Phys. Rev. 
B 81, 075124 (2010)] to the case of off-axis incidence. For the general case where both the refractive 
indices and layer thicknesses are random, we obtain the long-wave and short-wave asymptotics of 
the localization length over a wide range of incidence angles (including the Brewster "anomaly" 
angle). At the Brewster angle, we show that the long- wave localization length is proportional to the 
square of the wavelength, as for the case of normal incidence, but with a proportionality coefficient 
substantially larger than that for normal incidence. In mixed stacks with only refractive-index 
disorder, we demonstrate that p-polarized waves are strongly localized, while for s-polarization the 
localization is substantially suppressed, as in the case of normal incidence. In the case of only 
thickness disorder, we study also the transition from localization to derealization at the Brewster 
angle. 



PACS numbers: 42.25.Dd,42.25.Fx 



I. INTRODUCTION 

Anderson localization is one of the most fundamental 
and fascinating phenomena in the physics of disordered 
systems. Predicted in the seminal paper of Anderson 1 
for spin excitations, and extended to the case of elec- 
trons and other one-particle excitations in solids (see, 
for example, Refi^) and also applied to classical waves^, 
this very general phenomenon has become a paradigm of 
modern physics 4 .. Despite considerable efforts, the theo- 
retical framework of Anderson localization in higher di- 
mensions (D > 1) is far from complete-, especially in the 
case of classical waves where the effects of absorption 6 , 
gain 7 



' - and polarization^— are significant. 



In contrast, the one-dimensional case (D = 1) has 
been studied extensively for both quantum mechanical 
and classical waves (see, e.g., Refs^i^). In the systems 
with short-range correlated disorder, it is known that 
all states are localize d 13 ' 14 . One of the main manifesta- 
tions of localization is the exponential decay of the am- 
plitude of a wave propagating through an infinite disor- 
dered sample. This decay is the result of the interference 
of multiply scattered waves, and its spatial rate is called 
the Lyapunov exponent, 7, whose inverse value 7 _1 is 
a characteristic length describing localization in an in- 
finite sample. By itself, however, the reciprocal of the 
Lyapunov exponent does not provide comprehensive in- 
formation about the transport properties of disordered 



media for all cases (see Ref^ for details). Moreover, it 
is unlikely that this quantity can be measured directly, 
at least in the optical regime. 

A further manifestation of localization is the exponen- 
tial decay of the transmission coefficient of a long, finite 
sample. The characteristic length of this decay is the 
transmission length It , also denoted by I n in the case of 
a sample of N layers. In the localized regime, where this 
length is much smaller than the sample size, in general 
the transmission length coincides with the localization 
length I, and also with the inverse Lyapunov exponent 

The Anderson localization of classical waves in one- 
dimensional disordered systems has been studied in 
detail^^ 1 ^— ) and it has been shown that in the long 
wave region where the interference is weak, the localiza- 
tion length demonstrates a universal behavior, growing 
in proportion to the square of the wavelength, i.e., I oc A 2 . 

In recent years, we have witnessed the rapid emergence 
of a new field of research in metamaterials — artificial ma- 
terials which exhibit a negative refractive mdex-£~— . In 
such materials, the wave vector k, the electric field vec- 
tor E and the magnetic field vector H form a left-handed 
(LH) coordinate system, in contrast to the right-handed 
system (RH) that is applicable to normal or regular ma- 
terials; for this reason, metamaterials are sometimes re- 
ferred to as left-handed materials. In metamaterials, the 
directions of the phase velocity and the energy flow are 
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opposite. This feature can strongly affect Anderson lo- 
calization in metamaterials. Indeed, in stratified me- 
dia formed of alternating layers of normal materials and 
metamaterials, the phase accumulated during propaga- 
tion through a right-handed layer diminishes in its prop- 
agation in a left-handed layer— and, as a consequence, in- 
terference will be weakened and localization suppressed. 

Already, one of the first study22 of Anderson localiza- 
tion in the presence of metamaterials has revealed the 
striking behavior that, in the particular case of an alter- 
nating stack of right- and left-handed layers of the same 
thickness and randomly varying refractive indices, the lo- 
calization is strongly suppressed. Its functional form in 
the long- wave region changes from the standard behavior 
of I oc A 2 to I oc A 6 and, subsequently, it was shown 2 ^ that 
in such stacks the localization length differs from the in- 
verse of the Lyapunov exponent — the first such example 
of this surprising behavior. 

While the disorder is one-dimensional, the random 
stacks are actually three-dimensional objects and so the 
vector nature of the propagating field and, in particular, 
its polarization can strongly influence localization, lead- 
ing to its suppression or even complete derealization. 

For stacks comprising only normal layers, the effects of 
polarization have been extensively studied. The problem 
has been studied using an approach based on stochas- 
tic differential equations^, with excellent agreement ob- 
tained at short wavelengths (for normal incidence), and 
also in the long wavelength limit, between the theoreti- 
cally predicted localization length and that obtained by 
direct simulation. The derealization associated with the 
Brewster angle anomaly has been studied in the long 
wave limit using an effective medium approach^, and at 
short wavelengths in the framework of a random phase 
approximation^. Numerical simulations for the off- axis 
case are given in Refi 2 ^ while, in Refj 2 ^, the polarization 
properties of localization have been considered experi- 
mentally. 

Although the approaches described in RefaS r 11 ' 2 ^ 26 
give rich information about the effects of polarization on 
the properties of localization in normal materials, a gen- 
eral expression for the localization length, applicable for 
broad range of input parameters (including at the Brew- 
ster anomaly angle) is still missing. Furthermore, there 
are no results available for the effects of polarization on 
localization in the presence of metamaterials. 

In this paper, we extend our earlier study 2 ^ to the case 
of off-axis incidence and provide a comprehensive analyt- 
ical and numerical treatment of the effects of polarization 
on the localization length I. We consider both homoge- 
neous stacks, formed by only normal material layers or by 
only metamaterial layers, and also mixed stacks formed 
by an alternating sequence of right- and left-handed lay- 
ers with random thicknesses and refractive indices. We 
derive explicit asymptotic expressions for the localization 
length at the short and long wavelength limits that are 
applicable for the Brewster angle, and demonstrate their 
excellent agreement with direct simulation. 



The paper is organized as follows. In Sec. [HI we de- 
scribe the model and outline the theoretical treatment. 
The derivation of the asymptotic forms for the localiza- 
tion length at short and long wavelengths is presented in 
Sec. IIIII In Sec. IIV1 we present the results of numerical 
simulations for the localization length for s- and p- polar- 
izations in both homogeneous and mixed stacks. Here, we 
adopt the conventional definition for polarization, with s- 
and p-polarization referring respectively to the cases in 
which the electric and magnetic fields are perpendicular 
to the plane of incidence. Finally, in Sec. [V] the depen- 
dence of the localization length on the angle of incidence, 
at a fixed wavelength, is considered. 

II. THEORETICAL STUDIES 
A. Model 

We consider the transmission and localization proper- 
ties of a one-dimensional, multi-layered, disordered stack 
which consists of N layers composed of either right- 
handed or normal (r) materials, left-handed (I) meta- 
materials, or mixed stacks comprising alternating layers 
(r and I) of each (see Fig. []}. 
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FIG. 1: (Color on line) The geometry of the structure. 

All layers are statistically independent and, in the most 
general case, the thickness of each layer, and its dielectric 
and magnetic permittivities, are random quantities with 
given probability densities. All lengths in the problem 
are measured in the unit of the mean layer thickness and, 
therefore, are dimensionless quantities. 

The main subject of our interest is the transmission 
length^ 

N 
(ln|2V|) 

where is the transmission coefficient of the N layer 
stack for a plane wave with a given incidence angle (rela- 
tive to the surface normal of the stack) . Angular brackets 
are used to denote averaging over realizations of all ran- 
dom parameters. In the limit of a stack of infinite length, 
i.e., N — > oo, the transmission length coincides with the 
localization length, i.e., 

I = lim l N . (2) 

N-HX> 
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The calculation of the transmission length (fTJ) requires 
the transmission coefficient Tjv of the TV-layer stack for a 
wave of a given polarization and a given incidence angle. 
Such a calculation can be based on the transfer matrix 
methodic, the interface iteration method^, and the layer 
iteration method 2 ^. We choose the last of these and build 
on the treatment that was used successfully in our pre- 
vious study 24 for the case of normal incidence. 

The method is based on the exact iteration of the re- 
currence relations 



Tn 



T r > — 1 tr> 



1 - R n -ir n ' 

R f 2 
1 

1 - R n -ir n 



(3) 
(4) 



for the total transmission T n and total reflection R n co- 
efficients of a n- layer stack, n = 1, . . . , N, in which both 
the input and output media are free space, and with ini- 
tial conditions set to To = 1 and Rq = 0. Here t n and r n 
are the transmission and reflection coefficients of the n th 
layer, with layers being enumerated from n — 1 at the 
rear of the stack through to n = N at the front. 

Such an approach is quite general and is applicable to 
an arbitrary choice of polarization (s) or (p), the angle 
of incidence 9, the type of layer (r or I), the wavelength 
A, and the amount of absorption or gain. All these input 
data are incorporated into the transmission and reflection 
characteristics t n and r„ of a single layer. 

The recurrence relations ([3]) and (|4]), together with the 
definitions (fTJ) and © enable the numerical calculation 
of the transmission length In in the most general case. 
If the calculated transmission length / n is much smaller 
than the stack length N, and is independent of N, then 
In may be identified as the localization length, i.e., In — 
I. 

In the work reported in this paper, we deal with a spe- 
cific model in which the dimensionless thicknesses d of 
each layer are independent, identically distributed, ran- 
dom variables, d = 1 + 5^, where 5d is uniformly dis- 
tributed in the range [— Qd,Qd] with < Q4 < 1. The 
dielectric and magnetic permittivities of the layers are 
represented in the form 



e = ±(l + 5 u ) 2 , fi = ±l, 



(5) 



where the upper and lower signs respectively correspond 
to a normal material or a metamaterial. The random 
part 8 V of the refractive index, 



±(1 



(0) 



is uniformly distributed in the range [— Q V ,Q U ], < 
Q v < 1. Accordingly, the model takes into account both 
refractive index disorder and layer thickness disorder. 



B. Theoretical analysis 

Our theoretical treatment is based on a weak scatter- 
ing approximation (WSA) in which the magnitude of the 



reflection coefficient of each layer |r„| <C 1 is the primary 
small parameter in the theory. To understand when this 
condition is valid and may be used, we next consider ex- 
plicit expressions for the reflection r and transmission t 
coefficients of any single layer: 



t 



(W 2 )e 



i/3 



p 2 e 2iP ' 



(7) 



1 _ p2^P > I 

where p is Fresnel interface reflection coefficient given by 



P 



Z cos 9 U — cos 9 
Z cos 9 V + cos 9 ' 



Z 



z-\ 
z, 



s — polarization 
p — polarization. (8) 



In these equations, 



kdv cos 9,, , cos 9,. 



k = 2ir/\, (9) 



A is the dimensionless free space wavelength, and 



1 



> 



(10) 



is the layer impedance relative to the background (free 
space). 

There are two cases in which the magnitude of the sin- 
gle layer reflection coefficient is small, i.e., |r| -€i I. The 
first is characterized by weak refractive index disorder 
Q v <C 1 and corresponds to the incidence angle 9 being 
smaller than its critical value 



= shr 1 (I -Q v ). 



(11) 



In this case it is the small magnitude of the Fresnel re- 
flection coefficient \p\ -€i 1 which leads to a small, single 
layer reflection coefficient (|r| ^ 1) at all wavelengths. 
The second case corresponds to the long wavelength limit 
(A > 1) in which the single layer reflection coefficient 
is small (|r| <^ 1) for an arbitrary incidence angle 9 
due to the asymptotically small value of the multiplier 
1 1 — e 2l/3 1 oc <C 1 which appears in the expression for r 
in Eq. ©. 

Within the WSA approximation, we commence the 
derivation of general forms with the linearized recurrence 
relations ([!]) 

lnT„ = lnTi,„_i + lni„ + i?„_ir„, 
Rn = r n + R n ^t 2 nl (12) 

and solve them to yield 

N N N j-1 

\nT N = j2 lnt J + J2J2 r i~ m+ir i II *p' ( 13 ) 

j — 1 m—2j—m p=j — rn+2 

from which we may compute ensemble averages. 

In what follows, we summarize the key theoretical 
results 2 -^ applicable to mixed and homogeneous stacks. 
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Mixed stacks 



2. Homogeneous stacks 



A mixed stack, which hereafter is abbreviated by M- 
stack, is composed of alternating layers of right-handed 
(r) and left-handed (I) materials (see Fig.[IJ. Within the 
model, there are simple relations that may be derived^ 
between the averaged values of an analytic function g of 
the transmission and reflection characteristics of single 
layers of left- and right-handed materials, 



(g(t r )) = [g{ti)Y, (g(r r )) = {g(n)r 



(14) 



As a consequence, the transmission length of a M-stack 
depends only on the properties of a single right-handed 
layer and is expressible in terms of only three averaged 
quantities: (r) 2 , (lnt), and (t 2 ), in which the subscript r 
(referring to a right-handed layer) has been omitted: 



Here. 



1 _ 1 
In I 



-Re (lnt) 



b I J ^ \ lm. 



\(r) 2 \+R C ((r) 2 (t 2 y) 



is the inverse localization length, 



(15) 



(16) 



1 



2/l ri 



b I 1 — exp(— 2/l m ) 

/ |(r) 2 |+Re((r) 2 (t 2 )*) 

^ i-K* 2 >l 2 

is the inverse ballistic length, and 

= l-cxp(-z) ^ 



\(r) 2 



(17) 



(18) 



To characterize the transition from localization to bal- 
listic propagation, we introduce the crossover length for 
a M-stack2i 



1 



In |(t 2 } 



(19) 



and two characteristic wavelengths Ai and A2 defined by 



N = l(\ 1 (N)), N = l(X 2 {N)). 



(20) 



For long wavelengths, that part of the spectrum for which 
A <C Ai(AT) corresponds to the localization regime, for 
which In = I. In turn, the wavelength range A 3> A2(iV) 
corresponds to ballistic propagation and In — b. The 
region Ai < A < A2 is the transition region from localiza- 
tion to ballistic propagation. In what follows, we charac- 
terize these regions with long wavelength asymptotes for 
the localization length I and the crossover length I. 



The homogeneous stack (abbreviated from here as a 
H-stack) is composed of layers of the same type of mate- 
rial (either r or I). By virtue of the symmetry relations 
dH]), the results for statistically equivalent stacks of ei- 
ther normal materials or metamaterials are identical. 

The transmission length of a H-stack is then 



111 



(r) 



(l-(t 2 )) 2 _ 

while the inverse localization length is given by 



1 



-Re (lnt) - Rc 



(rf 



l-(t 2 )' 



(21) 



(22) 



It follows from Eq. (|21[) that the crossover length lh of 
the H-stack is 



l h = - 



1 



ln<t 2 )| : 



(23) 



characterizing the transition from the near ballistic 
regime to the far ballistic regime of propagation. 



III. ASYMPTOTIC BEHAVIOR OF THE 
TRANSMISSION AND LOCALIZATION 
LENGTHS 

The results (fTSJ) - fI7]). (|2Tj) . ([2"2")l for the transmis- 
sion, localization and ballistic lengths of homogeneous 
and mixed stacks are quite general. In this Section we 
apply them to the specific model described in Subsection 
III AI We emphasize that within the model both two types 
of disorder are present and are taken into account in all 
intermediate calculations. However, in all final results we 
keep only the leading terms. The higher order corrections 
with respect to the relative perturbations Q Ut d of the re- 
fractive index and thickness distributions are generally 
omitted. 



A. Short-wave asymptotics of the transmission 
length 

When the incidence angle is smaller than the critical 
angle 8 C (fTTj) , the short wavelength asymptotics of the 
localization length can be easily obtained from Eqs (fl())) 
and (f2"2"| . At short wavelengths, the phase of the field is 
a strongly fluctuating, random quantity so that (t 2 ) s» 0, 
(r) « (p), and (In « (ln(l — p 2 )). As a consequence, 
the localization length for both mixed and homogeneous 
stacks takes the form 



-Mi-p 2 ))-(p) 5 



(24) 
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For s-polarization, the logarithmic term in Eq. 
always dominates and so 



Ql 



12 cos 4 ' 



(25) 



while the second term in ([24]) provides a higher order 
correction of order O (Qi) ■ For p-polarization, however, 
the corresponding correction cannot be omitted since the 
first term vanishes at the Brewster angle — it/ A. Thus, 
for p-polarization, 



569 
~96~ 



1 

7 

43 



11 cos 2 (29) 
12 cos 4 9 



Qi 



120 cos 8 9 



8 cos 20 + — cos 40 
8 



cos 69 H cos i 

96 



.(26) 



Accordingly, at the Brewster angle, the localization 
length is given by 



/ = 



45 



(27) 



Note that in the case of normal incidence, the results (|2"5)) 
and (|26l) coincide with those presented in Refs&^ i 15 ' 16 . 

If the incidence angle 9 exceeds the critical angle 9 C 
(fTTj) . then total internal reflection occurs (i.e., the mag- 
nitude of the Fresnel reflection coefficient becomes unity) 
and so the WSA fails in the short wave region. If the in- 
cident angle is sufficiently far above the critical value C , 
then the exponent 2i/3 in Eq. {2} is real and negative 
and thus the magnitude of the single layer transmission 
coefficient is exponentially small. As a result, we obtain 
the following expression for the transmission length 

— = Im(/3) = fclm(d\/sm 2 0-^ 2 }, for sin0 > 1-Q V . 
In 

(28) 

The right hand side of this equation is independent of 
the stack length N and so it formally coincides with the 
inverse localization length I. However, its origin is related 
to attenuation by tunneling, rather than to Anderson 
localization. In the attenuation regime, the transmission 
length does not distinguish the left- and right-handed 
layers since it depends on the square of the refractive 
index v and so is the same for equivalent I or r layers. 
Moreover, it does not distinguish the polarization of the 
light. 

The average of Eq. (|2"5)) can be calculated in closed 
form for the uniform distribution of the refractive index 
and we obtain the short wave asymptotic of the transmis- 
sion length in the attenuation regime from the expression 



1 

In 



4Q„ 



sin 

2 



l-Q„ 

sin0 



- (i-<9.) 



(i-0,) 2 



(29) 



and see that it is proportional to the wavelength. We 
see also that it holds for both polarizations, and also for 
both H- and M-stacks. 



The main contribution to this length in the short wave 
region coincides with that of the first terms in Eqs (fTo]) 
and ([22]) for the localization lengths of H- and M-stacks. 
These terms in Eqs (ITB1) and (|2"2"|) dominate in the short 
wave region, as we will demonstrate below, and thus these 
terms give the correct values for the transmission length 
in the attenuation regime (for both short wave and long 
wave regions) despite the applicability of the WSA being 
violated. 



B. Long- wave asymptotics for homogeneous stacks 

1. Homogeneous stacks: s-polarization 

For long wavelengths, the transmission length can be 
deduced from the general result in Eq. (|2"2"1) . In this limit, 
the mean values of (ln£), (t 2 ) and (r) that enter (|2"2"|) for 
s-polarization take the form 



(ln£) = ik cos 9 + 



ikQl k 2 Ql 



6 cos 9 6 cos 2 9 



Q 2 , 

1 + T 



k 2 Qt 



(r) 



1 + 

40 cos 2 0^3 
k Q v 

6 cos 9 6 cos 2 9 



(30) 



(2- 



cos 2 9 



3Qj 



x (1 + p 
ik^Ql 



(4- 



9 cos 
2 cos 20 
2<:os 3 I 3 + 



cos 



10 ' 

[i + Q 



15 



(t 2 ) 



*k 3 Ql 
56 cos 3 

1 + 2ifccos0 

2fc 2 Q 2 



(1 + Qd)> 



(31) 



1 



■ cos 



] ikQl 
3 cos 
k 2 Ql 
3 cos 2 
3k 2 Qj 
20 cos 2 



- 2fc 2 cos 2 
(3 + 2 cos 2 ( 



1 



(32) 



Applying the results of Eqs (|3"01) - ([52"|) in the expression 
for the transmission length (|2T|) . we obtain 



1 

Ln 



k 2 Ql 
6 cos 2 

Qi 

72 cos 4 



Ql 



15 3 cos 2 



sin 2 (fcA cos0) 
N 



(33) 



which is correct to the order of Q^. This expression can 
be further simplified given that Q u <C 1 and Qd <1. In 
this approximation, the transmission length takes form 



1 



IN 



k 2 Ql 
6 cos 2 



NQl /sin(/WVcos0) 



12 



kN cos 



(34) 



6 



from which it follows that the localization length is given 
by 



I = 



3A 2 cos 2 f 



I < N. 



(35) 



From this, it is seen that the localization length of a ho- 
mogeneous stack in the long wavelength limit does not 
depend on the thickness disorder for s polarization. We 
verified this through exact numerical calculations (sec 
FigHJ. 

Eq. ([35[) reproduces the result first derived in Ref. 9 
and, for normal incidence, coincides with the localization 
length derived more recently i n 18 ' 24 . For wavelengths 
such that I > N > A/ cos 0, there is a near ballistic regime 
where the ballistic length b is the same as the localization 
length 



b = 



3A 2 cos 2 



I > N > 



A 



cost 



(36) 



The crossover length in the far ballistic region (|23|) is 

A 



h = 



Air cos 9 ' 

while the ballistic length bf in this region is 



3A 2 cos 2 , 
2tt 2 Q 2 



1 



NQl 
12 



From Eqs (J35j>, ([37l) and ([20]), it follows that 

2N 
3~' 



Ai 



COS 

A 2 = A-k Ncos 6, 



(37) 



(38) 



(39) 
(40) 



for all reasonable values of the input parameters such that 
Ai < A2. For A < Ai, waves are localized with the local- 
ization length given by Eq. ([3"5[) . The wavelength range 
Ai < A < A2 corresponds to the near ballistic regime in 
which the ballistic length b is given by Eq. (f36|) . The 
longer wavelengths, A > A2, correspond to the far ballis- 
tic regime in which the ballistic length bf is given by Eq. 



2. Homogeneous stacks: p -polarization 

In the case of p— polarized waves, the mean values of 
(lnt), (t 2 ) and (r) 2 , at longer wavelengths, take the form 



(lnt) = ifccos0 + 



ifcQ 2 cos30 ifcQ 4 sin 2 ' 
6 cos 2 9 2 cos 9 



-(i + Ql) 



k 2 Qj cos 2 29 f Q 2 
6 cos 2 9 V 3 

k2Q l (2-3cos20) 2 + ^ 
4Ocos 2 v ; V 3 



k 2 Ql tan 2 9 



(2-3cos6»), 



(41) 



(r) 



ikQl(2 - cos 29) ikQism 2 9 



(1 + Ql) 



6 cos 9 2 cos 9 

^%2(3 + 1Ocos20 - cos 40) (l + % 



24 cos 2 9 
k 2 Qt 



,f31- 52 cos 20 + 17 cos 40) f l + % 
80 cos^ 9 \ 3 



fc 2 Q^tan 2 



(2 - 3 cos 20) +ak 3 



(42) 



(r } = 1 + 2ik cos 9 - 2k 2 cos 



2 2 a , ^Q? c °s30 



6 cos 2 



k 2 Ql 
3 cos 2 

cos 2 I 



-(2cos 2 + 3cos 2 20) 



(43) 



where the expression for the coefficient of the cubic term 
a in ([421 ) is given in the Appendix. The transmission 
length In, given by Eq. ([2"T[) . can then be expressed in 
the asymptotic expansion form as 



In 



k 2 Q 2 v cos 2 29 f 1 t Q 



_ n 

6 cos 2 



24 cos 4 



+ 



+ 



k 



2nd 



17280 cos 6 

Qt_ 

187V cos 4 V2 



sin 2 ( k N cos 0), (44) 



where 



19 7 19 

hAO) = 1 cos 20+ —cos 40+ —cos 60, (45) 

6 15 30 

h 2 {9) = 1415 - 1664cos20 - 188cos40 

+ 512 cos 60+ 141 cos 80. (46) 

This expression can be further simplified given Q v <C 1 
and Qd <C 1. We obtain the final form of the transmission 
length correct to the order of 0(Q 4 ): 



J_ = fc 2 g 2 cos 2 20 / Q\ 
In 6cos 2 V 3 



k 2 Qt 
24 cos 4 



cos 2 sin 2 (fcA^cos0),(47) 



187V cos 4 V2 
with the localization length being given by 

1 _ fc 2 Q 2 cos 2 20 k 2 QlQ 2 cos 2 29 



I 



+ 



6 cos 2 
k 2 Qj 
24 cos 4 



18 cos 2 



hi(9), 



(48) 



where h\{9) is given by (1451) . We have verified numeri- 
cally the asymptotic formula (j4"B"[) in Sec IV AI and found 
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that it is in excellent agreement with the exact numerical 
calculations for angles of incidence of up to 80°. 

Equation (j48j) generalizes the corresponding expression 
for the localization length obtained in Ref.— and is ap- 
plicable to incidence at the Brewster anomaly angle of 
9 = 7r/4. At this angle, the first two terms in Eq. 
vanish, leading to 



I = 



45A 2 
l6iT 2 Qt 



7T 

4' 



(49) 



This length is proportional to Q~ 4 , in contrast to the Q~ 2 
dependence applicable for any incidence angle away from 
the Brewster angle. However, its wavelength dependence 
of order O (A 2 ) remains the same for all angles less than 
the critical angle. In the case of weak disorder (Q„ <C 
1), this means that the localization length can be made 
arbitrarily large at the Brewster anomaly angle relative 
to that realizable at other incidence angles. 

Then, using Eq. (|43|) in Eq. (|23|) . we may deduce 
that the characteristic wavelength A2 is identical to that 
obtained for s-polarization (|4"0)) . Similarly, Ai, obtained 
using Eq. (|48l) . is given by 



Ai 



cos 9 



cos 2 29 



cos 2 9 



hx(6). 



(50) 



The significance of the threshold wavelengths Ai and A2 
is the same as for s-polarization. Waves with A < Ai 
are localized, while the wavelength range Ai < A < A 2 
corresponds to the near ballistic regime. Similarly, the 
ballistic length b has the same form as the localization 
length (|49| 



45A 2 
16tt 2 Q 4 



7T 

4' 



(51) 



Thus, the transition from localization to the near ballistic 
regime takes place without any change of the scale. The 
wavelength range A > A 2 is the far ballistic region in 
which the ballistic length 6/ is given by 



1 

bj 



k 2 Ql cos 2 29 I Ql 
6 cos 2 9 V 3 
Nk 2 Qj (3 _ 2 
lScos^U C0S 



(52) 



We emphasize that the results obtained in this subsection 
are applicable only to homogeneous stacks composed of 
normal material or metamaterial layers. 



C. Long-wave asymptotics for mixed stacks 

1. Mixed stacks: s-polarization 

Substituting the asymptotic forms ([3T))) - (l3"2"T) into Eqs 
TBI and (|17[) . we derive an expression for the reciprocal 



transmission length: 

1_ k 2 Q 2 , 
In 



1 1 - f(Na s ) 



where 



3cos 2 6> \2 3 + Ccos 4 i 
k 2 (l 2 

- ^-(3 + Ccos 4 0), 



3 cos 2 9 ' 

the function / is as defined in Eq. flTBl . and 



2QI 
Ql 



(53) 



(54) 



(55) 



Eq. (|53[) describes the transition from localization to 
ballistic propagation at long wavelengths and, in the limit 
as N — > 00, we obtain the following expression for the 
localization length 



/ 



3A 2 cos 2 6> 3 + Ccos 4 



2vr 2 Q 2 1 + Ccos 4 



(56) 



The ballistic length formally corresponds to the opposite 
extreme, i.e., as N — > 0, 



b = 



3A 2 cos 2 9 
2tt 2 Q 2 ' 



(57) 



and coincides with the result for a H-stack in s- 
polarization. 

The characteristic wavelengths Ai and A2 take the form 



Ai 
A 2 




2N /l + Ccos 4 <9 
3V 3 + Ccos 4 6»' 



cost' V 3 



(58) 
(59) 



Again, for the range A < Xi(N), waves are localized, 
while ballistic propagation occurs for very long wave- 
lengths A > A 2 (^V). The transition wavelengths Ai j2 are 
of the same order and the intermediate region Ai (iV) < 
A < A 2 (iV) corresponds to the crossover between local- 
ization and ballistic propagation. 



2. Mixed stacks: p -polarization 

While we have derived a general expression for the 
transmission length that is applicable at arbitrary an- 
gles of incidence for a M-stack in p-polarization, its form 
is quite complex and so it is presented only in the Ap- 
pendix. In what follows, we look at a number of partic- 
ular cases. 

For incidence at angles away from the Brewster angle, 
the transmission length, according to Eq. (IA.1I) . is given 
by: 



1 

In 



k 2 Q 2 v cos 2 29 (\ 
3 cos 2 9 V 2 



1 - f(Na p ) 



2 + cos 2 26» + C cos 4 9 / 



(60) 
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where 



k 2 Q? 
3 cos 2 



-(2 + cos 2 20 + Ccos 4 



(61) 



The localization length may be deduced from Eq. ([60 
by taking the limit as N — s> oo, i.e., 



3A 2 cos 2 9 2 + cos 2 29 + ( cos 4 9 
2tt 2 Q 2 cos 2 29 cos 2 29 + < cos 4 6 ' 



(62) 



Correspondingly, the ballistic length may be obtained 
by calculating the limit as N — > in Eq. ([BTJl) : 



6 = 



3A 2 cos 2 
2tt 2 Q 2 cos 2 29' 



(63) 



The transmission length for the Brewster anomaly an- 
gle can be deduced by substituting 9 = 7r/4 in (|A.1[) : 



1 

In 



4k 2 Qj 
45 



V 



121Q 2 



44 



60 



-/(^« P ) 



(64) 



/ 



Here, the second term in parentheses is always smaller 
than the first and therefore, at the Brewster angle, the 
transmission length as a function of the wavelength ex- 
hibits the same dependence, i.e., 



In = I = b = 



45A 2 



(65) 



and is independent of the length of the stack. As in 
the case of s-polarization, the transmission length in the 
localized regime for p-polarization behaves as O (Q^ 4 ) 
and exceeds the localization length far from the Brewster 
angle. Note that the localization length (|65p for the M- 
stack at 9 — ir/4 is the same as for the homogeneous 
stack (|49|). 

The transition between the localized and ballistic 
regimes is again described by two characteristic wave- 
lengths Ai and A2: 



Ai 



A 2 




cos 2 29 + C cos 4 9 
cos 2 26> + Ccos 4 




a/4 + 2 cos 2 29 + C cos 4 9. 



j, (66) 



(67) 



The expression for Ai (1661) is obtained under the assump- 
tion that the angle of incidence is sufficiently far from 
the Brewster angle. At the Brewster angle, Ai and A2 
are given by 



Ai 



4ttQ 2 



N 



A 2 = 4nQ L 



16' 



(68) 



(69) 



Note that A2 3> Ai. This means that at the Brewster 
angle, waves such that A < Ai(iV) are localized, while 
the ballistic region Xi(N) < A becomes divided into two 
sub-regions. The near ballistic sub-region is bounded by 
the two characteristic lengths, i.e., Xi(N) < A < \2(N), 
and the far ballistic region corresponds to very long waves 
^(N) < A. As we explained previously, the localization 
length and the ballistic lengths in each of the two ballistic 
subregions are described by the same expression (|65l) . 



IV. RESULTS OF NUMERICAL SIMULATIONS 

We now present results of our comprehensive numeri- 
cal study of the properties of the transmission length as 
a function of wavelength and angle of incidence. Results 
are presented for direct simulations based on exact recur- 
rence relations ©-©j and are compared with those ob- 
tained from the analytic forms ([I5 ]) -(jT7 ) . (p2"Tj) . and (f2"2"1) . 
and their short and long wavelength asymptotic forms 
derived in Sec. IIIII 



A. Homogeneous stacks 

1. Subcritical angle of incidence 

We consider transmission through a H-stack charac- 
terized by the parameters: Q v = 0.1, Qd = 0.2 at the 
incidence angle 9 = 45°, which is less than the critical 
angle 9 C — sin _1 (0.9) » 64.16° and coincides with the 
Brewster anomaly angle for a layer with a mean refrac- 
tive index of v = 1 . 

We begin with s-polarization, and consider a stack of 
length N — 10 4 , using N r = 10 4 realizations for ensemble 
averaging. For the given parameters, the characteristic 
wavelengths are Ai ~ 36 and A2 ~ 8.9 x 10 4 . Plotted in 
Fig. [2] is the transmission length as a function of wave- 
length. Fig. Ufa) corresponds to relatively short wave- 
lengths A < 10 2 and represents mainly the localized part 
of the spectrum where N > Ijy « /. Fig. [2^b) corre- 
sponds to longer waves and mostly displays the ballistic 
part of the spectrum where In > N. In both panels, 
the red solid lines represent Zjv(A), obtained by exact nu- 
merical simulation, and the blue dashed lines display the 
analytical form ([2~Tj) . The excellent agreement between 
these two curves for all wavelengths (in both panels) is 
evident. 

The curves displayed in Fig. [2f a) explicitly confirm the 
coincidence of the long wave asymptotes of the transmis- 
sion length in both the localized (A < 36, In — I) and 
near ballistic (36 < A, In — b) regions. The slanted, 
dashed line corresponds to the asymptotic forms ([33]) and 
([31>| while the horizontal dashed line corresponds to the 
short wave asymptote (|25|) . The corresponding curves of 
Fig. [UJb) display the transmission length in the ballistic 
regime. The near ballistic region, where In — b, occurs 
for 36 < A < 8.9 x 10 4 , while the transition to the far 
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FIG. 2: (Color online) Transmission length In versus wave- 
length A for Q v — 0.1, Qd = 0.2 and 9 = 45° for s-polarized 
waves; panels (a) localized part of the spectrum and (b) bal- 
listic part of the spectrum. Red solid curve: numerical simu- 
lation; Blue dash curve analytic form (|21l) . 



FIG. 3: (Color online) Transmission length In versus A for 
Q v = 0.1, Qd = 0.2 for p-polarized waves at the Brewster 
anomaly angle 8 = 45 . Panel (a) localized part of the spec- 
trum; panel (b) the transition from localization to ballistic 
propagation. Red solid curve: numerical simulation; blue 
dash curve: analytic form (|21l) . 



ballistic region, where In = bf, occurs for A w 8.9 x 10 4 . 
The upper and lower dashed lines respectively display 
the near and far ballistic lengths of Eqs (j5oT) and ([3"5|). 
We observe that the results for s-polarization are entirely 
consistent with those reported previously for the case of 
normal incidence 2 ^. 

Figure[3]presents the corresponding results for the case 
of p -polarized waves. Here, we consider a much longer 
stack of N — 10 6 layers, the characteristic wavelengths of 
which are Ai « 19 and A 2 ~ 8.9 x 10 6 . Fig. displays 
the transmission length spectrum for comparatively short 
wavelengths A < 10 2 and corresponds mainly to the 
localized part of the spectrum where N > In w I. The 
results of Fig. [3Jb) correspond to longer waves A > 10 2 
and display the ballistic part of the spectrum A > Ai 
where In > N. In both panels, the red solid and the 
blue dashed lines respectively display Zjv(A) obtained by 
exact numerical calculation and the analytic form (1211) . 
and we see that their agreement is excellent. 

The results of Fig.[3Ja) show that the long wave asymp- 
tote of the transmission length in both the localized re- 
gion A < 19 (where In — I), and the near ballistic region 
19 < A < 10 2 (where In = b) coincide exactly. The 
slanted dashed line corresponds to the asymptotic form 



(|49| while the horizontal dashed line represents the short 
wave asymptote . We observe that the short and long 
wave limits for the localization length at the Brewster 
angle are proportional to Q^ 4 and hence are two orders 
larger than in the case of s-polarization. These numerical 
results confirm the analytical results presented earlier in 
Sec. lIHB2l 

Fig. |3]Jb) characterizes the ballistic regime which com- 
prises a near ballistic region (19 < A < 8.9 x 10 6 ) in which 
In = b, and a far ballistic region where In — bf, with the 
transition between the two occurring at A w 8.9 x 10 6 . 

The upper and lower dashed lines respectively display 
the asymptotes for the near (p)Tj) and far (|52[) ballistic 
lengths. The ballistic length, over the entire ballistic re- 
gion, including the transition from the near to far bal- 
listic regime, is very well described by Eq. (|47l) . The 
oscillatory nature of this transition is due to Fabry- Perot 
resonances between the first and the last interfaces of the 
stack and is much more pronounced than for the case of 
s-polarization. We observe that the envelope of the trans- 
mission length curve is confined from below by Eq. (|47[) 
in which the sine term is replaced by unity. This is the 
long dashed black curve of Fig. Elb) . While this highly 
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oscillatory region also occurs for s-polarized waves, it is 
not apparent in Fig. [5] since the chosen stack length 
(N = 10 4 ) was not sufficiently long to exhibit the fea- 
ture. 

In the samples with only refractive index disorder (i.e., 
Qd = 0), the localization length displays strong oscilla- 
tions at intermediate wavelengths 0.3 < A < 2 for both 
polarizations. Equation (|2"2"|) is also in an excellent agree- 
ment with the numerical calculations. We also note that 
the thickness disorder smears out these oscillations, with 
only few oscillations remaining for Qd = 0.2 (see Figs 

2. Supercritical angle of incidence 

When the angle of incidence exceeds the critical angle, 
i.e., 8 > 8 C = sin _1 (l — Q v ), the exponential wave decay 
can be attributed not only to Anderson localization but 
also to attenuation inside the individual layers. 
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FIG. 4: (Color online) Transmission length In of the homo- 
geneous stack with N = 10 4 , Q v = 0.1, Qd = 0.2 versus 
wavelength A for s-polarized waves at the supercritical inci- 
dence angle 9 = 75°. Red solid curve: numerical simulation; 
Blue dash curve analytic form (|21|) . 

In Fig. |H we plot the transmission length spectrum 
for a s -polarized wave in which the parameters of the 
problem are the same as for Fig. [21 apart from the angle of 
incidence which is 8 = 75°. In this case, the characteristic 
wavelengths are Ai as 99 and A2 as 3.2 x 10 4 . Since, as 
noted previously in Sec. IIII A[ Eq. (|2"Tj) can serve as a 
good interpolation formula for the transmission length 
in the short and long wave regions, we have plotted the 
results of the exact numerical simulation (red solid curve) 
together with those predicted by Eq. (|2"Tj) (blue dashed 
curve) to demonstrate the quality of the agreement. 

In Fig. SJ the dotted line displays the short wavelength 
asymptotic (f29|) . the black dashed, slanted line displays 
the localization length (|35p . and the dashed dotted line 
shows the far ballistic asymptotic given by (|38p . At short 
wavelengths, A < 2, the form of the transmission length 



spectrum is determined mainly by attenuation or "tun- 
neling" effects. The transmission length is proportional 
to the wavelength and is well described asymptotically 
by Eq. the dotted black slanted line in Fig. g) 

For longer waves, the form of the transmission length is 
the same as is observed below the critical angle, and is 
described well by Eq. (|3~4"]) . 

Anderson localization is realized only in the interme- 
diate wavelength region, 2 < A < 99, with In as I with 
the localization length given by Eq. (1351) . and shown in 
the black dashed slanted line of Fig. 

For p-polarization, the spectral behavior of the trans- 
mission length is qualitatively equivalent to that for s- 
polarization and so we do not present this here. There is 
excellent agreement between the exact numerical calcu- 
lation, the theoretical result of Eq. (|2"2"|) . and the short 
([2~4l and long (l4"5)l wave asymptotic forms. 



B. Mixed stacks 

1. Subcritical angle of incidence 

We first consider the case of s-polarized wave propa- 
gation through a mixed (i.e., alternating layers of normal 
and meta-materials) stack of length N = 10 4 . The pa- 
rameters are the same as those adopted in Sec. IIV A 11 
i.e., Q v — 0.1, Qd = 0.2, N r = 10 4 , and the incidence 
angle is 6 = 45°, which is less than the critical angle 
8 C = sin _1 (0.9) = 64.16° and coincides with the Brew- 
ster angle for the single layer with mean refractive index 
v = ±1. For these parameters, the characteristic wave- 
lengths given by Eqs (155)) and ([59")) are Ai as 28, and 
A2 as 115 correspondingly. 

In Fig. [5j the red solid line and the blue dashed line re- 
spectively display results from the numerical simulation 
and the analytic form (based on the WSA) for the trans- 
mission length as a function of wavelength, with these 
two coinciding to high accuracy. 

The form of the transmission length spectrum is sim- 
ilar to that observed for the case of normal incidence^. 
The short wavelength asymptotic form shown as 

the horizontal dashed line in Fig. [5j is the same as for a 
H-stack. For A < Ai = 28, all waves are localized, with 
the localization length (156")) shown by the upper dashed, 
slanted straight line. The transition from localization 
to ballistic propagation occurs in the wavelength range 
Ai < A < A2 and is well described by Eq. ([551 . Ballistic 
propagation occurs for A > A2 = 115 and is character- 
ized by the ballistic length (|57l) which differs from the 
M-stack localization length ([53]) . in contrast to the case 
of H-stacks. 

Figure [5] displays the transmission length spectrum 
for a M-stack of length N = 10 6 in p-polarized light, 
with all other parameters identical to that for the s- 
polarization simulations. In this case, the characteristic 
wavelengths are Ai as 19 (from Eq. ([6"6")) ) and A2 ~ 12007T 
(from Eq.([67|). 
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FIG. 5: (Color online) Transmission length versus A for a 
M-stack in s-polarized light with Q v — 0.1, Qd = 0.2 and 
N = 10 4 for 6 = 45°. Red solid curve: numerical simulations; 
Blue dash curve: analytic form (|15|) . 
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FIG. 6: (Color online) Transmission length versus A for a 
M- stack in p-polarized light with Q v — 0.1, Qd = 0.2 and 
N = 10 6 , at the Brewster angle 8 = 45° red solid line. The 
blue dashed line shows results for s-polarization and a H- 
stack, re-plotted for comparison. 



2. Supercritical angle of incidence 
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FIG. 7: (Color online) Transmission length versus A for a Mi- 
stack in s-polarized light with Q v = 0.1, Qd = 0.2 and N = 
10 4 , and for the supercritical incidence angle 9 — 75°. Red 
solid curve: numerical simulations; Blue dash curve: analytic 
form ([Toll . 

We now consider a case in which the angle of incidence 
6 = 75° exceeds the critical angle 6 C = sin (1 — Q v ) = 
64.16° (jlip. In Fig. [3 we present the transmission length 
spectrum for s-polarized light and display results from 
the exact numerical calculation (red solid line) and the 
analytic form (long dashed blue curve). The expectation 
that Eq. (fT5|) would serve as a good interpolation formula 
for the transmission length in the short and long wave re- 
gions, as anticipated in Sec. IIII Al is borne out by the 
results of Fig. [7] The short wave (dashed dotted line) 
asymptotic form (|29|) and the long wave (black dashed 
line) asymptotic form (|57[) respectively coincide with the 
numerical results for A < 1 and 200 < A. In the inter- 
mediate region 1 < A < 200, however, the theoretical de- 
scription underestimates the actual transmission length 
since the WSA is no longer valid for the chosen, super- 
critical angle of incidence. For p-polarization, the results 
are qualitatively the same, but with the discrepancy at 
the intermediate wavelengths even more pronounced. 



The results of the numerical simulation and the WSA 
analytical forms ([T5|) , (fTfl)) coincide and are displayed by a 
single red solid line. Localization occurs for A < Ai « 19, 
while the transition from localization to ballistic prop- 
agation occurs at A ~ Ai. In contrast to the case of 
s-polarization, the transition is not accompanied by a 
change of scale and is given by the same wavelength de- 
pendence (|6"5|) . The same asymptotic (|65l) also holds for 
wavelengths A > A2 ~ 12007T, which defines the transi- 
tion from the near to the ballistic regime. As a conse- 
quence of the disorder Q u = 0.1, the short wave local- 
ization length (|24l) (horizontal dashed line) is two orders 
of magnitude larger than that for s-polarized light (cf. 
Fig©. 



C. Mixed stacks with refractive-index disorder 

In our earlier paper—, we demonstrated that at nor- 
mal incidence a disordered mixed stack, with only refrac- 
tive index disorder, could substantially suppress Ander- 
son localization. Indeed, the suppression is so strong that 
even the usual quadratic dependence on wavelength (i.e., 
0(A 2 )) of the localization length at long wavelengths was 
shown to change to 0(A 6 ). In contrast, the introduction 
of the thickness disorder in combination with the refrac- 
tive index disorder induces strong localization at long 
wavelengths, with the localization length returning to its 
expected quadratic dependence on wavelength^ 4 .. In this 
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section, we consider the effects of polarization on long 
wavelength localization in M-stacks. 
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FIG. 8: (Color online) Transmission length In versus A for a 
M-stack with Q v = 0.25, Qd = and 8 — 30° for p-polarized 
light (cyan dashed dotted curve, N = I0 6 ) and s-polarized 
light (red solid curve, N — fO 5 ; green dashed curve, N = 10 7 ; 
blue dotted curve, N = 8 x 10 8 ). 



Figure [5] displays transmission length spectra for a 
mixed stack with only refractive index disorder for an an- 
gle of incidence of 6 = 30°. Four curves are displayed: for 
p-polarized light and a stack of length N = 10 6 (dashed 
doted cyan curve), and for s-polarized light and three 
stacks of lengths N — 10 5 (solid red curve), N = 10 7 
(dashed green curve) and N = 8 x 10 s (blue curve). 
There is a striking difference between the two polariza- 
tions: in the case of p-polarized light, there is strong 
localization at long wavelengths (A < 10 2 ), with the lo- 
calization length showing 0{\ 2 ) dependence; in contrast, 
the localization length for s-polarized light is much larger 
and shows the 0(X e ) dependence as occurs for normal in- 
cidence. Note that for s-polarization, the localization re- 
gions in Fig. [5]are bounded from above by the wavelength 
limits A < 5,9, and 12 for stacks of length N = 10 5 , 10 7 , 
and 8 x 10 8 respectively. 

This asymmetry between the polarizations suggests 
that the suppression of localization is due not only to 
the suppression of the phase accumulation but also to 
the vector nature of the electromagnetic wave. Because 
of the symmetry of Maxwell's equations between the elec- 
tric and magnetic fields, it is to be expected that for a 
model in which there is disorder in the magnetic perme- 
ability (with e = ±1) the situation will be inverted, with 
localization for p-polarized waves being suppressed and 
with s-polarization showing strong localization. 

In concluding this section, we emphasize that the del- 
icate phenomenon of the suppression of localization oc- 
curs only for refractive index disorder, and that the in- 
troduction of any thickness disorder leads to the strong 
localization (see Sec. IIVB~T1 and Ref^). 



V. TRANSMISSION LENGTH AS A FUNCTION 
OF THE INCIDENCE ANGLE 



A. Homogeneous stacks 
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FIG. 9: (Color online) Transmission length In versus in- 
cidence angle 6 for a homogeneous stack with Q v — 0.1, 
Qd = 0.2 for (a) A = 0.1 (upper panel), (b) A = 10 (lower 
panel). Red solid curve: numerical simulations; Blue dash 
curve analytic form (|22[l . The top set of curves in each of the 
panels are for p-polarization while the bottom set of curves 
are for s-polarization. 

We next consider the angular dependence of the trans- 
mission length of a homogeneous stack for a given wave- 
length. As in earlier simulations, we work with the pa- 
rameters Q v = 0.1, Q d = 0.2, and N = 10 6 . Figure M 
displays the transmission length as a function of the an- 
gle of incidence 9 for both s- and p-polarizations. In 
each panel (upper: A = 0.1, lower: A = 10), the solid 
red curve displays the results of the numerical simula- 
tion while the blue dashed line corresponds to the WSA 
analytic form (|2"2"|) . with the top and bottom sets being 
for p- and s-polarizations respectively. 

For the short wavelength A = 0.1 (Fig. (Ha)), 
the analytic form agrees perfectly with the simulations. 
While for s-polarization, the transmission length de- 
creases monotonically with the angle of incidence, the 
transmission length for p-polarization displays a pro- 
nounced maximum at the Brewster anomaly angle (at 
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9 sa 46° for these parameters). In the supercritical 
regime, 9 > 9 C « 64°, attenuation is the dominant mech- 
anism for localization and hence the behavior of the two 
polarizations coincide. 

For long wavelengths, as in Fig. HJb), we see that for 
extreme angles of incidence (e.g., for 9 > 80° for the 
wavelength A = 10), the theoretical prediction departs 
markedly from the simulation results. Similar departures 
for intermediate wavelengths (e.g, for A = 1) also exist 
for angles of incidence 9 > 85°. 

We have also calculated the localization length for the 
very long wavelength of A = 40 as a function of the angle 
of incidence (for the same parameters as for FigJS]), for 
which the expansion (|4"8f is applicable. There is excellent 
agreement between the exact numerical calculation and 
the asymptotic form (|48p for angles of up to 80°. (Since 
this plot is very similar to FigJ^b), it is not included in 
the manuscript.) 



B. Mixed stacks 
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FIG. 10: (Color online) Transmission length In versus inci- 
dence angle 6 for a mixed stack with Q v = 0.1, Qd = 0.2, 
for (a) A = 0.1 (upper panel), and (b) A = 1 (lower panel). 
The top and bottom curves are respectively for p- and s- 
polarizations. 

We now consider the angular dependence of the trans- 
mission length for mixed stacks and, in Fig. llOl we plot In 



as a function of the angle 9 for a stack of length N = 10 6 
at the two wavelengths A = 0.1 (Fig. ITUT a)) and A = 1 
(Fig. [TCT b)). In either case, the calculated transmis- 
sion length does not exceed the stack length and so, for 
subcritical angles, our calculations display the true lo- 
calization length. For the shorter wavelength A = 0.1, 
the form of the transmission length for both polariza- 
tions is similar to that observed for homogeneous stacks 
(cf. Fig. 02a)), and we also note that the analytical form 
(JTHJ) agrees perfectly with the results from the numerical 
simulations. 

Fig. HOf b) displays results for an intermediate wave- 
length A = 1 with the lower solid red and blue dashed 
curves respectively displaying the results of numeri- 
cal simulations and analytical predictions (1161) for s- 
polarization, (bottom curves), while the upper solid green 
and brown dashed curves display simulations and analyt- 
ical predictions (|16|) for p-polarization. The agreement 
between simulations and the theoretical form is again ex- 
cellent for angles of incidence less then the critical angle, 
9 < 9 C , while for angles greater then the critical angle, 
the discrepancies that are evident are again explicable 
by the breaking down of the WSA at extreme angles of 
incidence. 



C. Alternating homogeneous stacks 

In this section, we present an example of true delocal- 
ization arising from the vector nature of the electromag- 
netic field. This was first pointed out by Sipe et aJiS, 
in which an analysis applicable at long wavelengths was 
presented. More recently, the analysis has been extended 
to short wavelengths^. The condition for the Brewster 
anomaly can be satisfied for a homogeneous stack (i.e., 
with all layers being either normal materials or all being 
metamaterials) with only thickness disorder, and with 
alternating refractive indices (i.e., with refractive indices 
va and vb respectively for odd and even numbered lay- 
ers). 

We proceed in a similar manner to that of RefJ£ and 
consider a stack in vacuum with va = 1 and vb = 1.5, 
and with layers whose random thicknesses are uniformly 
distributed in the interval d G [0.8,1.2] (i.c, Q d = 0.2). 
We note that in the case of p-polarization, the applica- 
bility of the weak scattering approximation is heightened 
in the vicinity of the Brewster angle since each layer is 
almost transparent. 

In Figs [11] (a) and (b), we respectively plot the trans- 
mission length as a function of the angle of incidence at 
a short wavelength A = 0.1, and also at an intermedi- 
ate wavelength A = 0.5. The lower and upper curves 
are respectively for s- and p-polarizations, and we see, 
somewhat surprisingly, that the numerical simulations 
and the analytic forms obtained within a WSA frame- 
work are essentially identical for both polarizations for 
arbitrary incidence angles, apart from the discrepancies 
evident for extreme angles 9 > 87° at A = 0.5. The sur- 
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FIG. 11: (Color online) Derealization at the Brewster angle 
6b', (a) at a short wavelength A = 0.1, and (b) an intermediate 
wavelength A = 0.5. Top and bottom curves are respectively 
for p- and the s-polarizations. 



prising element is that the theoretical framework based 
on the WSA appears to work over a much wider range 
of angles and polarizations than that suggested by strict 
validity of the WSA. In the figure, the theoretical de- 
scription for s-polarization (green curve) overlays the re- 
sults of the numerical calculation (red curve). The same 
is true for p-polarization, with the theoretical prediction 
(black dotted line) overlaying results from the numerical 
calculation (cyan solid curve). 

In these calculations, the stack length was N = 10 4 and 
so waves are delocalized for incidence angles 55° < 6 < 
59° around the Brewster angle Ob = arctan(1.5) « 56.19° 
where the transmission length In > N. The localiza- 
tion properties for the corresponding homogeneous stack 
composed of metamaterial layers with vb = —1.5 is the 
same as that shown in Figs [TT] for normal layers with 
vb = 1.5. The WSA based theory also appears to work 
over a reasonably broad range of wavelengths, although 
for the intermediate wavelengths 10 < A < 50 there are 
some differences between simulations and the theoretical 
prediction. 



VI. CONCLUSIONS 



We have investigated the effect of polarization on 
the Anderson localization in one-dimensional disordered 
stacks composed entirely of either right- or left-handed 
layers, as well as mixed stacks with alternating sequence 
of normal and metamaterial layers. 

Our analysis has generalized the results obtained 
earlier 24 for the case of normal incidence to the case of an 
arbitrary angle of incidence, with a particular attention 
paid to the localization at the Brewster angle. Based on 
this approach, we have carried out a comprehensive study 
of the localization length as a function of both the angle 
of incidence and the polarization of the incident wave for 
various types of disorder. 

In the case of general disorder, where both refractive 
index and thickness of the layers are random, we have 
derived the long- and short-wave asymptotics for the lo- 
calization length for a wide range of incidence angles, 
including the Brewster angle. At the Brewster angle, we 
have shown that the localization length continues to ex- 
hibit a quadratic dependence on wavelength (as in the 
case of the normal incidence), but that the coefficient 
of proportionality becomes paramctrically larger, being 
proportional to Q~ 4 , rather than Q~ 2 (Q <C 1), as for 
the case of the normal incidence. 

Our theoretical study not only characterizes the local- 
ization and ballistic propagation, but also describes per- 
fectly the crossover between these two regimes. We have 
also shown that the transition from localization to ballis- 
tic propagation in the vicinity of the Brewster angle in a 
mixed stack is given by a single scale (|65|) . In the case of 
thickness disorder, we have shown that, at the Brewster 
angle, Anderson localization is suppressed completely. 
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Appendix: Transmission length for p-polarization 
and mixed stack 



In this Appendix, we provide an asymptotic form for 
the transmission length in the case of p-polarization for 
a mixed stack. The reciprocal transmission length (|15[) 



15 



takes the form 



and 



1 

In 



+ 



+ 



where 
fi -- 

h -- 



7^ cos 2 20 fl + f 
6 cos^ \ 3 

-fi%-(2-3cos20) 2 
40 cos 2 

k 2 ^ 



4 ^(2-3cos20) (1 + ^ )tan 2 ' 



1+ fl 

Qd 



k 2 Ql 
48 

2fc 2 Q 4 cos 2 20 



3 cos : 



(43 - 55 cos 20) tan 2 

+ h + h + h + h + h 



2Q 2 V (2 + cos 2 29) + 4Q 2 cos 4 9 

2 Qt cos 2 20 

3co S 2 + ff 1 + 92 + 93 + ff4 + .95 



4Q 2 (2 + cos 2 20) +8Q 2 cos 4 



(A.l) 



k 2 QtQ 2 (61 25 



8 cos 2 
5 



19 



d ' ^ + ^ cos 20+ ^- cos 40 
8 6 3 



6 



cos 69 



cost 



k 2 Ql 
12 cos 2 V 120 
49 



24 

1927 121 



(A.2) 



cos 29 + cos 46* 
5 30 



— cost 
15 



24 



(A.3) 



k 2 QlQ 2 , /7013 1763 „„ 3772 
32^1 (w-^5- COS20+ T35- COS40 



1391 163 

cos b9 H cos < 

135 270 



(A.4) 



fc 2 Q» /633 19058 9521 

cos26»H cos46» 

32 cos 2 9 V 4 75 75 

^ cos 6^ + — cos 86) , (A.5) 
25 300 / ' y ' 

k 2 QtQi (3 + io C os20-cos40) 2 , (A.6) 



9i 



<)2 



93 



<Ji 



95 



2k 2 Ql 
45 cos 2 9 
11 cos 66») , 

k 2 QtQ 2 d /-.Q 

18cos 2 fl l 
cos 69) , 

k 2 Ql 
270 cos 2 9 

2362 fifl 
cos 69 - 



(41 - 77 cos 29 + 41 cos 40 



5 cos 20+ 10 cos 40 



1963 - 
121 



16534 



cos 20 



8968 



5 5 
k 2 QlQ 2 (247 
60 cos 2 V 3 



139 cos 20 



928 



■ cos 40 



103 



cos 60 



11 



|g| ( 3 + lOcos20-cos40) 2 . 



(A.7) 

(A.8) 
cos 40 

(A.9) 

(A.10) 
(A.ll) 



The expansion (|A.1[) is valid for any angle of incidence 
for mixed stacks at long wavelengths. The factor / is 
given by (|18p in this expression and characterizes the 
transition from localization to ballistic propagation. 

For the sake of completeness, we also provide in this 
Appendix the expression for the cubic coefficient a in the 
expansion (|12"]). 
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All expansions in (|30 [l - (j32|l . (|4T ]> -([43 j) and in the Ap- 
pendix can be readily obtained by using symbolic manip- 
ulation package such as Mathematica. 
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